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Abstract 

The spin-charge gauge approach to the metal-insulator crossover (MIC) and 
other anomalous transport properties in high-T c cuprates is briefly reviewed. A 
U (1) field gauging the global charge symmetry and an SU (2) field gauging the global 
spin-rotational symmetry are introduced to study the two-dimensional t — J model 
in the limit t > J. The MIC, as a clue to the understanding of the 'pseudogap' (PG) 
phase, is attributed to the competition between the short-range antiferromagnetic 
order and dissipative motion of charge carriers coupled to the slave-particle gauge 
field. The composite particle formed by binding the charge carrier (holon) and spin 
excitation (spinon) via the slave particle gauge field exhibits a number of peculiar 
properties, and the calculated results are in good agreement with experimental data 
for both PG and 'strange metal' phases. Connections to other gauge field approaches 
in studying the strong correlation problem are also briefly outlined. 



PACS numbers: 71.10.Pm, 71.10.Hf, 11.15.-q, 71.27.+a 



1 Introduction 



The discovery of high temperature superconductivity 20 years ago (Bednorz and Miiller 
1986) has posed a great challenge to the condensed matter physics community It is true 
that there is no consensus yet on theoretical interpretation of the normal state properties 
and superconducting mechanism for these cuprates. However, enormous progress has been 
made in material sample preparations, experimental probes of various physical properties 
using the finest tools and theoretical studies from many different points of view. Although 
physicists working in this field still have different opinions, most researchers would more 
or less accept the following scenario: 

These cuprates share a layered structure incorporating one or more copper-oxygen 
planes. The parent compounds (e.g. La 2 Cu0 4 ) contain one electron per site which should 
be metals according to the band theory, but are indeed Mott insulators due to strong 
Coulomb repulsion. Upon doping beyond a certain concentration (a few per cent) these 
compounds become superconducting. Meanwhile, in the temperature range above the 
superconducting transition (usually called pseudogap (PG) 'phase' on the underdoped 
side and strange metal (SM) 'phase' on the optimal, or slightly overdoped, side), various 
physical properties, including optical, magnetic and transport, are anomalous and cannot 
be described by the standard Landau Fermi liquid (FL) theory. Most people believe that 
these unusual properties are due to correlation effects in doped Mott insulators. One has 
to first study these properties using simple models, like the Hubbard model, or t — J 
model in the strong coupling limit, to see whether the new physics in the PG and SM 
phases beyond the FL theory can be essentially understood. Surely, there are other effects, 
like inhomogeneities and electron-phonon interactions, but most likely they are not the 
dominating factors contributing to the new physics. 

The strong-correlation approach was pioneered by Anderson (1987), reviving his ear- 
lier work on frustrated spin 1/2 system at triangular lattices (Anderson 1973, Fazekas 
and Anderson 1874). Instead of antiferromagnetic (AF) long range order (LRO), a spin 
liquid ground state consisting of singlet resonant valence bonds (RVB) was considered. 
Using the Gutzwiller-projected BCS-wave function Anderson argued that this RVB state 
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is naturally evolving into superconducting state upon doping. Later, the spin-charge sep- 
aration concept (Baskaran et al 1987, Kivelson et al 1987, Zou and Anderson 1988) was 
introduced, with a spinon carrying the spin degree of freedom and a holon carrying the 
charge degree of freedom. Using the slave boson mean field theory (MFT) (Kotliar and 
Liu 1988) and Gutzwiller approximation (Zhang et al 1988) a d-wave superconductivity 
and pseudogap state were predicted using the RVB approach. The gauge field fluctua- 
tions implementing the single-occupancy constraints beyond the slave boson MFT were 
considered by a number of authors (Baskaran and Anderson 1988, Ioffe and Larkin 1989, 
Nagaosa and Lee 1990). Very recently, a comprehensive review on the strong-correlation 
approach to high T c (HTS) cuprates, focusing on the gauge theory approach was provided 
by Lee et al (2006). 

For the last few years we have developed a different gauge approach to attack the HTS 
problem, using a (7(1) field to gauge the global charge symmetry and an SU(2) field to 
gauge the global spin symmetry ( Marchetti et al 1998, Marchetti et al 2000, Marchetti et 
al 2001, Marchetti et al 2004a,b, Marchetti et al 2005). Unlike the slave-boson approach, 
the charge degree of freedom in our approach is carried by a fermion, while the spin degree 
of freedom is represented by a boson. We found this approach to be very efficient in 
treating the most pronounced phenomenon in the PG phase, namely the metal-insulator 
crossover (MIC) as well as various transport, optical and magnetic properties. In this 
short review we will outline the basic ideas of our approach and summarize our main 
results in comparison with experiments. We will attempt to make the presentation more 
readable and emphasize the physical interpretation, referring readers to original papers 
for technical details. 

2 Metal-insulator crossover as a clue to the HTS 
problem 

The PG phase was anticipated by the RVB approach and has been studied very thoroughly 
in an enormous amount of experimental and theoretical work which has been very well 
summarized in recent reviews (Timusk and Statt 1999, Norman and Pepin 2003, Lee et 
al 2006). In our view the most spectacular phenomenon of the PG phase is the MIC 
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observed in underdoped, non-superconducting cuprates in the absence of magnetic field 
and a similar phenomenon in superconducting samples when an applied strong magnetic 
field suppresses the superconductivity. 

First of all, this MIC is a rather universal phenomenon. A minimum in resistivity 
(around 50 - 100K) and a crossover from metallic conductivity > at high temper- 
atures to insulating behavior ^ < at low temperatures has been observed in heavily 
underdoped LSCO (Takagi et al 1992, Keimer et al 1992, Ando et al 2001, Ando et al 
2004), non-superconducting Bi 2+a; Sr2_j / Cu0 6 ± < 5 (Fiory et at), non-superconducting YBCO 
(Wuyts et al 1996, Trappeniers et al 1999, Ando et al 1999, Ando et al 2004) and La-doped 
Bi-2201 (Ono and Ando 2003). It has also been observed in electron-underdoped NCCO 
(Onose et al 2001) and PCCO ( Fournier et al 1998). 

Moreover, such a MIC has also been observed in a number of superconducting samples 
when a strong magnetic field suppresses the superconductivity, including LSCO (Ando et 
al 1995, Ando et al 1996a, Boebinger et al 1996), La-doped Bi-2201 (Ando et al 1996a,b, 
Ono et al 2000), electron-doped PCCO (Fournier et al 1998) as well as in Zn-doped YBCO 
(Segawa and Ando 1999). 

This phenomenon has not had much attention so far, because many people attributed it 
to localization due to disorder. However, that interpretation is at odds with the following 
facts: (1) Including higher doping samples exhibiting MIC in strong magnetic fields, the 
estimate for kpl, where hp is the Fermi momentum and i the mean free path, ranges from 
0.1 to 25 at the MIC, i.e. from far below to far above the Ioffe-Regel limit {kpi ~ 1), 
characterizing the MIC due to disorder localization. (2) In LSCO samples with a — 
b in-plane anisotropy the MIC temperature of p a is different from that of pb (Dumm 
et al 2003) contradicting the 'unique' localization temperature, characteristic of the (at 
least standard) theory of localization. (3)A universality of suitably normalized resistivity 
(Wuyts et al 1996, Konstantinovic et al 2000) has been observed in terms of T/T*, where 
T* is roughly proportional to T at the MIC (Tmic) and can be identified as the PG 
temperature. All the above features are very difficult, if not impossible to explain from 
the localization viewpoint. On the other hand, as discussed below, the spin-charge gauge 
approach provides a rather natural interpretation for these unusual features (Marchetti 
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et al 2004a,b). 

Furthermore, once the MIC is recognized as an intrinsic property of cuprates, it is very 
difficult in conventional theory to reconcile the insulating behavior at low temperatures 
with the presence of a finite Fermi surface (FS), as shown by angle-resolved photoemission 
spectroscopy (ARPES) (See, e.g., Damascelli et al 2003). Meanwhile, these two seemingly 
contradictory phenomena can be easily accommodated in a slave-particle gauge theory, 
due to the Ioffe-Larkin (1989) rule, stating that the inverse conductivity of the electron 
is the sum of the inverse conductivity of the spinon and the inverse conductivity of the 
holon. This non-standard feature can be intuitively understood as a consequence of the 
gauge string binding spinon to holon: the velocity of the electron is determined by the 
slowest (not the fastest!) among spinon and holon. Then, if the fermionic excitation, 
spinon or holon depending on the approach, has a FS and therefore a metallic resistivity 
vanishing at T ~ 0, the electron can still have a metallic/insulating behavior at low T if 
the bosonic excitation does. If the leading contribution comes from the boson, without 
a detailed reference to FS, this can qualitatively explain the universality quoted in point 
(3) above. 

In our view the MIC is the clue to the understanding of the PG phase. The MIC in 
underdoped cuprates in the absence of magnetic field and MIC in superconducting samples 
when a strong magnetic field suppresses superconductivity is the same phenomenon with 
the same origin: as an outcome of competition between the AF short range order (SRO) 
and the dissipative motion of the charge carriers. We start from the Mott insulating state 
showing AF LRO. Upon doping beyond certain threshold the AF LRO is destroyed, being 
replaced by SRO, characterized by an AF correlation length £. Since holes distort the 
AF background and their average distance ~ intuitively, £ ~ where 5 is the 

doping concentration. This has been confirmed by the neutron scattering experiments 
(Keimer et al 1992). Our theoretical treatment proves £ (5\ ln5|) _1 / 2 , providing the 
first length scale, while the corresponding energy scale is the spin excitation (spinon) gap 
m s = J(S\ In <5| ) 1 / 2 , where J is the AF exchange interaction. A competing factor is the 
diffusive motion of the charge carriers with characteristic energy ~ Tm,h, where T is the 
temperature, while m h ~ 5 ft is the effective mass of the charge carrier (holon) in the PG 
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phase, t is the hopping integral. The corresponding length scale is the thermal de Broglie 
wave length \ T ~ (T5/t) -1 / 2 . At low temperatures £ < X T , the AF SRO dominates 
and the charge carriers become weakly localized (not exponentially) showing insulating 
behavior. We would like to emphasize that this 'peculiar localization' is mainly due to 
interaction rather than disorder, and it is also different from the standard Mott insulator 
which comes entirely from the Coulomb repulsion. On the contrary, at high temperatures 
£ > At, the diffusive motion of charge carriers prevails, exhibiting metallic conductivity. 
Therefore the competition of the real part of the spinon 'self-energy', the mass gap, and 
the imaginary part, the dissipation, gives rise to this spectacular phenomenon: MIC. 

It turns out that a number of experimental observations in the PG phase can also be 
explained by the competition of these two factors and the 'composite' nature of low-energy 
excitations. Roughly speaking, the spin (spinons) and charge (holons) excitations in the 
PG phase behave like 'separate particles' in their scattering against gauge fluctuations, 
which renormalizes their properties and dominates the in-plane transport phenomena. 
However, at small energy-momentum scale the gauge field binds spinon and anti-spinon 
into magnon 'resonance'. Similarly spinon and holon are bound into electron 'resonance' 
with non-Fermi-liquid properties, showing up in ARPES experiments as a 'quasiparticle 
peak' and responsible for interlayer transport. In particular, the 'composite' nature of par- 
ticles is exhibited in the anomalously large spectral weight away from the quasi-particle 
peak seen in ARPES (See, e.g., Nagaosa and Lee 1990, Lee and Nagaosa 1992, Laughlin 
1997, Orgad et al 2001). In a sense, the low-energy excitations in these strongly corre- 
lated materials are no longer well-defined quasi-particles of FL theory, but rather these 
loosely bound 'composite particles', taken as a first step going beyond the well established 
condensed matter physics paradigm. 
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3 Spin-charge decomposition: statistics and symme- 
tries 



Our basic assumption is that the low-energy physics of the Cu-0 layers in HTS cuprates 
can be described qualitatively by a two-dimensional t — J model: 

H = Pg H (-* Z) c ^' + J ^ ' ^ - \ n ^)) P Gi (3-1) 

<ij> a 

where i correspond to Cu-sites, a = 1, 2 is the spin index and Pq denotes the Gutzwiller 
projection, eliminating the double occupancy. Numerical simulations yield t ~ 0,4eV , 
J ~ 0.13el/ which we use in calculations of physical quantities. The qualitative behavior 
of the low-energy physics in the c-direction is obtained by adding an interlayer hopping 
term. 

Once the t — J model is accepted, a key question is how to handle the Gutzwiller 
projection. A straightforward way is either to use numerical techniques (Gros 1988, 1989, 
Paramekanti et al 2001, 2003), or MFT ( Zhang et al 1988, Anderson et al 2004). Another 
way is to decompose the electron into holon and spinon: 

c a = hs a (3.2) 

A simple counting of degrees of freedom (dof) shows that the decomposition (13.21) has 
an unphysical redundancy, i.e., a local gauge symmetry. In fact, the electron has 4 dof 
(c a , c*), while the spinon s a being still of two components, also has 4 dof, and the holon 
h being spinless has 2 dof. One has to impose the constraint of no double occupancy, 
which in slave particle formalism becomes holonomic. Then the redundant dof (4+2-1=5) 
corresponds to a local phase factor by which the holon can be multiplied, while the spinon 
is multiplied by its inverse: a U(l) gauge symmetry. It is very useful to make manifest the 
'hidden' gauge symmetry by introducing an emergent U(l) gauge field, coupling spinons 
to holons. 

3.1 Different choices of spin-charge decomposition 

Accepting the spin-charge decomposition as a good technique to explore the low-energy 
physics of cuprates, one still has several options for its implementation. The first choice 
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regards the statistics of spinon and holon. If no approximations are made all these choices 
are equivalent and a proof of this equivalence for some of them is given in Frohlich and 
Marchetti (1992). However, as soon as approximations are made, like mean field (MF), 
they become intrinsically different. Standard choices are: s a fermion and h hard-core 
boson or s a hard core boson and h fermion. These choices are at the basis of traditional 
slave boson (Zou and Anderson 1988, Nagaosa and Lee 1990, Lee and Nagaosa 1992) and 
slave fermion (Arovas and Auerbach 1988, Yoshioka 1989, Dorey and Mavromatos 1991) 
approaches. To describe PG the more sophisticated approach of slave boson pursued in 
Wen and Lee (1996), Lee et al (1998) is obtained by replacing a single boson h with a pair 
of bosons h a: enhancing the slave-particle symmetry from U(l) to SU (2). A slave- fermion 
version of this is discussed in Farakos and Mavromatos (1998). 

Moreover, in 2D there are also other 'anyon' statistics available, in particular the 
'semion statistics' advocated originally by Laughlin (1988), in which an interchange of the 
fields produces a factor e ±l 5. The product of two semions is still a fermion: (e ±l ^) 2 = —1. 
Notice also that for semions a kind of 'hard-core exclusion' holds (Frohlich and Marchetti 
1992), here ensured by Pq- 

To adequately describe the FS of cuprates determined by ARPES (Damascelli et al 
2003) is an important guideline in making choices: for optimally doped materials one 
finds a FS agreeing with band calculations. This is clearly compatible with s a being 
a fermion, as in the slave boson approach, but it is troublesome for the slave fermion 
approach because h being spinless the MF Fermi momentum is expected to be doubled 
w.r.t. to the electron case (Lee et al 2006). For a semion, although the situation is not 
completely elucidated (Wu 1994, Haldane 1994, Chen and Ng 1995), it appears that a 
kind of generalized Pauli-principle holds, related to Haldane statistics (Haldane 1991), 
stating that for low T one can accommodate at most 2 (spinless) semions in the same 
momentum state and the semion distribution function is approximately twice the fermion 
distribution function. Hence a gas of spinless semions of finite density should give a 'FS' 
coinciding with that of spin \ fermions of the same density. We take here this generalized 
Pauli principle for granted. 
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3.2 Hint from the ID t — J model 

The ID t — J model in the limit J/t \ has been solved exactly by the Bethe Ansatz 
method (Ogata and Shiba 1990), while the correlation functions have been calculated 
using the conformal field theory (CFT). Applying the spin-charge decomposition approach 
we have reproduced the main results of these techniques within a kind of MFT (Marchetti 
et al 1996). A coordinate representation of the solution exhibits the following features: 

(1) Spin-charge separation: the charge and the spin dof are characterized by different 
physical behavior, in particular their velocities are different v s ^ v c (Luttinger liquid) 

(2) The low-energy physics of the charge dof is described by a free spinless fermion (holon), 
while the low-energy physics of the spin dof is described by an spin | AF Heisenberg model 
in a squeezed spin chain, obtained by omitting the unoccupied sites of the original ID 
lattice, where the spins are given by S = S^^-Sp with S a a 'renormalized' gapless spinon. 
S a can be identified as a Gutzwiller projected fermion in the squeezed chain (Marchetti 
et al 1996). The scaling limit of such model is the 0(3) nonlinear cr-model with vacua 
at = 7r, the value when gapless excitations appear in ID. 

(3) The electron field is a product of a holon and a spinon field together with a non-local 
'string phase': c ax = h x exj)[^Hj <x hjhj]Sax- It can be shown (Marchetti et al 1996) that 
both S ax and 

h. r <' : --- ■ !r . h - (3.3) 

describe semions. One can prove that S a is just a MF approximation of a non-local spinon 
field with a 'string phase' analogous to that in (13.31) . constructed out of bosonic spinons 
s a , and its explicit expression is given in Marchetti et al (1996). The semion statistics is 
crucial for obtaining the correct correlation functions in the scaling limit (Marchetti et al 
1996) (see also Ha and Haldane 1994). 

One should note that this statistics is compatible with an optimization of energy: in 
the semion approach the spinon s a appears in the hopping term in the Affleck-Marston 
(AM 1988) form AM=^ a s* ai s a j while in the AF-Heisenberg term in the RVB form, 
RVB=^ Q/3 e a [3S ai spj, as |RVB| 2 . This occurs because for a holon-empty site j the spinon 
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contribution turns out to be given by 








(3.4) 



where the spin flip induced by a x is due to an emergent AF structure. However, for 
a holon-occupied site there is an additional spin flip and this is the case for the final 
site of a hopping link. The peculiar feature of this spin flip occurs only in a SU(2)- 
formalism. The interesting point is that the following identity holds for s a hard-core 
boson: |AM| 2 +|RVB| 2 =1. Optimization of the t term suggests | <AM> |=1 but due to 
the previous identity this choice optimizes also the positive J term, since it corresponds 
to | <RVB> |=0. 

From (I3.3P one can see that a semion can be constructed out of fermions (and in 
similar way also from hard-core bosons) through a generalization of the Jordan- Wigner 
transformation (Fradkin 1991) with an exponent which is half of the standard one. The 
continuum analogue of (13.31) is 



where H(x) is the Heaviside step function. The representation (13. 5p makes it clear that 
the 'statistics transmutation' (from fermion to semion) in ID is obtained by adding an 
exponential of a 'step function' or 'kink' average of the density. 

3.3 Generalization to two dimensions 

If we accept the suggestion coming from the ID model for the statistics of holon and 
spinon, one should search for a semionic representation of the electron field in 2D. A 
semion field in 2D is constructed out of a fermion as in (13.51) by replacing the ID 'kink' 
H(x — y) with a 2D 'vortex' arg (x) = arctan(^) in the average of the density. 
It is convenient to rewrite the resulting exponential in terms of a gauge field 




(3.5) 




(3.6) 



as exp[i B(y) ■ dy\ 

To show that indeed B = {B^ fi = 1,2} is the gauge field associated with a vortex 
let us compute its field strength, using V x V arg (x — y) — S(x — y) we get 
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where p is the density of h (Notice that from V • V arg (x — y) = it follows V • B = 0, 
i.e. B is in the Coulomb gauge). Equation (j3.7p is the same as that appearing in the 
fractional quantum Hall effect (FQHE) (see, e.g. Fradkin 1991) replacing ~ with ~ where 
v is the filling at the plateau. Hence the vortices appearing in (13. 6p are analogous to those 
introduced by Laughlin in the FQHE and in fact a semionic representation of the electron 
was advocated by him (Laughlin 1988) in the early days of HTS. 

Actually the physics here is dual to that of FQHE: there the flux was integer and 
the charge fractional, while here it is the opposite. The above ideas have been made 
precise, within the spin-charge gauge approach, in Marchetti et al (1998), where the 
following decomposition of the electron has been proved, in terms of a fermionic holon 
and a bosonic spinon s a satisfying the Gutzwiller constraint ^ a s* s a = 1: 

c x ~ h^Sfvydffp^nvW)^^ (3. 8 ) 

where B and V are U(l) and SU{2) gauge fields, respectively, and P is a path-ordering 
needed for the non-abelian nature of V^. In ( 13.81) we have 

B{V) = ± £ Vaxgtf - j)(l - h^iy )), (3.9) 
j 

3 

V(y) = £> a V ar 9(y~m - ^^(2/°))(^i<^ai(2/ )) (3.10) 

a=l 

with s a j = s a j if j is on the even Neel sublattice and s a j = e a ps%^ if j is on the odd 
Neel sublattice (compare with (13.41) ). Here an AF structure emerges formally and the 
basic related assumption is that fields s a have a good continuum limit, but not s a . The 
physical spin turns out to be 

Sj ~ ~s* aj a aP s Pj = (-l) b \s* aj a a(3 s Pjl (3.11) 

where |j|=0 (1) for a site in the even (odd) Neel sublattice, so that the AF ordering 
naturally evolves as 5 \ to the LRO of the 2D AF Heisenberg model, in agreement 
with experiments in undoped cuprates. In (I3.8P both h x e % $* B ' d ^ and (Pe 1 ^ Vd ^) a i3Sj3 X 
are semion fields. 
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3.4 Our choice of mean field approximation 

Where do the U(l) and the SU(2) groups of (13.91) . ( 13.1 Oft came from? They are just 
the global charge- U(l) and spin-SU(2) symmetry (not to be confused with the SU(2) 
slave-particle symmetry of (Wen and Lee 1996, Lee et al 1998)) of the t — J model. In 
Marchetti et al (1998) it is shown that if these spin-charge symmetries are made local by 
introducing a coupling with a £7(1) field B and a SU(2) field V given by ( 13. 91) . then the 
obtained 'spin-charge gauged t — J model' is strictly equivalent to the original t — J model 
S. Of course, as soon as MF approximations are made, this approach becomes distinct 
from e.g. the slave boson or slave fermion, becoming a new 'slave-semion' approach. A 
variant of this is obtained by replacing SU(2) of spin with the U(l) subgroup unbroken 
in the AF phase. This is the choice made by Weng and collaborators (Weng et al 2000, 
Weng 2003). 

However, in the 'semionic' approach it is too hard to keep charge, spin and slave 
particle symmetries all exact and an approximation has been made to perform explicit 
computations. In the spin-charge gauge approach pursued by us, the spin-charge gauge 
symmetry is treated in a sort of MFT while the slave particle symmetry is kept exact. 
An opposite choice has been adopted by Weng (Weng et al 2000, Weng 2003), i.e. the 
charge and the (reduced) spin symmetries are kept exact as much as possible, but the 
slave-particle symmetry is treated in MFT. 

The second major choice in the MF approaches is the spinon order parameter: if AM- 
like i.e., < s* ai spj > (Affleck and Marston 1988) and nonvanishing, it leaves unbroken 
the global slave-particle symmetry but breaks the global spin SU(2); if RVB -like i.e. 
< So/3 e af3S a iS/3j > (Wen and Lee 1996, Lee et al 1998) and non-vanishing, it breaks the 
global slave particle symmetry but preserves the global spin SU(2). 

In the standard slave-boson approach to PG one has <RVB>^ 0, though in the 
S£/(2)-slave boson approach a U(l) subgroup is unbroken. A kind of slave-boson approach 
where in PG one has<AM>^ is pursued in Feng et al (2004). 

The choice of the MF order parameter has relevant physical consequences, because the 
Ioffe and Larkin (1989) rule does not appear naturally compatible with a MF in which 
path integral formalism this is obtained within a Chern-Simons theory 
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the slave-particle symmetry is broken. 

Since in the spin-charge gauge approach we have <AM>^ 0, but <RVB> =0 in MF in 
both PG and SM 'phases' the Ioffe-Larkin rule holds. Furthermore since in MF | <AM> 
| = 1, the semion approach seems to optimize both the t and the J terms, as discussed in 
ID (but here a rigorous proof is lacking) in compatibility with an AF structure. 

4 Spin-charge gauge approach: the key ingredients 

To summarize, the key physical ingredients of the spin-charge gauge approach of MFT 
to characterize the 'normal state' within the t — J model are: the Gutzwiller projection, 
tackled with spin-charge decomposition, the two-dimensionality needed for semion statis- 
tics, finally leading to a spinon gap, and the bipartite lattice structure needed for both 
AF and the ir — flux crossover, which is also a peculiar 2D phenomenon. 

4.1 Energy optimization and n — flux crossover 

Here we focus on the n — flux crossover. A key step is a 'semiclassical' treatment of the 

— * 

interaction of B in the holon hopping term th*hjU < ij > , where f/<jj>is a complex abelian 
gauge field given by 

U<H> = e^^^sUPe^^^Us^. (4.1) 

If p denotes a plaquette and < ij >G dp the links in the boundary of p, the argument of 
n<ij><=dpU<ij> plays the role of a magnetic flux. 

It has been rigorously proved by Lieb (1994) that at half-filling (5 = 0) the optimal 
configuration for a magnetic field on a square lattice in 2D has a flux n per plaquette 
at arbitrary temperature. On the other hand, it is well known that at low densities 
and high temperatures the optimal configuration has flux per plaquette. At T = 
it has been proven that the ground state energy has a minimum corresponding to one 
flux quantum per spinless fermion (Bellisard and Rammal 1990). Numerical simulations 
(Qin 2000) suggest that increasing T gives rise to a competition between these minima 
and at sufficiently high T only and tt flux survive. We expect that the perturbation 
introduced by the J term in the t — J model changes the boundary, but not the essence 
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of the phenomenon of the tt — > crossover itself, while further enhancing the 'disorder' 
lowers the temperature where only tt and flux states survive. 

In the spin-charge gauge MFT it is assumed that the region of tt flux of U corresponds 
to PG and that of 0-flux to SM, thus the 'melting of the 7r-fiux lattice' underlines the 
physics of crossover between the two 'phases', which in particular produces a modification 
in the structure of the FsJ 2 ) 

The optimization considerations discussed above suggest a MF approximation for B 
given by 

B(y) — B MF {y) = ± ^ V arg(y - j), (4.2) 

3 

which gives flux tt per plaquette. Furthermore, to obtain the tt — ► crossover in the MF 
treatment of the holon hopping one imposes, by a suitable spinon redefinition, that 



<ij>£dp 



1 PG 

e in SM 



so that the effect of Bmf on holons in SM is cancelled by the spinon contribution. The 
unbalanced presence of Bmf in PG produces the Hofstadter phenomenon: it converts the 
spinless holon h of SM with dispersion 

uj ~ 2t[(cos k x + cos k y ) — 5] (4.3) 

and Fermi momenta fcjr ~ 1 — 8, i.e. a 'large FS' roughly consistent with band calcula- 
tion, into a pair of 'Dirac fields', with pseudospin index corresponding to the two Neel 
sublattices and dispersion relation: 



uj ~ 2t(- v /cos 2 k x + cos 2 k y - 5) (4.4) 

restricted to the magnetic Brillouin zone, thus yielding a 'small FS' centered at (±f, if) 
with Fermi momenta kp ~ 5. As we will see, many crossover phenomena from PG to SM 
in transport physics are due to this change of holon Fermi momenta. 

2 Let us remark that even in principle not all phenomena usually referred to as 'pseudogap' should fit 
in this interpretation, but only those referring to the T* temperature discussed below. 
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4.2 Generation of the spinon mass 

Comparing (14.21) and ( 13. 9ft one can see that the MF approximation for B consists in 
neglecting fluctuations of h; it is therefore reasonable in a MF treatment of V to neglect 
the spinon fluctuations i.e. 

V(y) — > K F {y) = 5 a3 J2(-l) ]jl Vwg(y-j)(l - h*hj). (4.5) 

j 

The first term of V^ 1F can be gauged away , while the second term describes spin vortices, 
centered at the holon positions, of opposite vorticity (or chirality) for holons on two Neel 
sublattices. These chiral spin vortices are reminiscent of those introduced by Shraiman 
and Siggia (1988) in analyzing the single hole problem, but here their density is finite 
5 > 0. 

Due to their AF nature the spinons in this approach have a low-energy effective action 
described by a 0(3) non-linear a-model plus an interaction with spin vortices appearing 
in Vmf- Due to vortex chirality, the spatial average of the interaction term linear in 
Vmf vanishes, hence in MF approximation the leading spinon- vortices interaction term is 
given at large scales by < V^ IF ■ Vj^p > s* a s a where <> denotes spatial average []. Naively 
one expects this average to be proportional to 5, as this is the density of holons, hence 
of vortices. However, due to the long-range tail of vortices one actually gets ( Marchetti 
et al 1998) a logarithmic correction, which, we will see later on, is physically relevant. 
Therefore, as a result of the interaction with vortices, in MF the spinons acquire a mass 
gap m s given by 

m 2 s ~ < vl IF .vl IF> ~5\\ Yl 5\. (4.6) 

The gap generation here is a kind of wave localization phenomenon: The gapless waves 
(here spinons) interacting with a finite (supercritical) density of impurities (here the 
vortices), are getting localized, acquiring a gap. The AF massive structure of the spinon 
in our approach prevents boson condensation at low T, in contrast to the slave-boson 
approach. 

3 Sclf-consistently the AF LRO is assumed absent 
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Since the spinon gap has an obvious relation with the AF correlation length, £af, 
its appearance implies that the finite (supercritical) density of vortices convert the long- 
range AF of the undoped (5 = 0) material into a short-range AF and the derived doping 
dependence (roughly £mf ~ 5~ l l 2 ) is consistent with the experimental data in neutron 
experiment in LSCO (Keimer et al 1992). Actually, by including the spinon thermal 
mass m T of the renormalized classical treatment of the 0(3)-model both the order of 
magnitude and the qualitative 5 — T dependence of the neutron data (Keimer et al 1992) 
are reproduced. The T independent term found in experiments which is hard to explain in 
terms of a standard 0(3) model, as in the spin-fermion approach (Abanov et al 2003), is 
here exactly due to the introduction of the spin vortices, a unique feature of this approach. 

4.3 Reizer singularity of the slave particle gauge field 

As already remarked, in slave-particle approaches a slave-particle U(l) gauge field couples 
spinon and holon. In the spin-charge gauge approach this can be identified with the gauge 
field, Ap, of the 0(3)-model rewritten in the CP 1 from, i.e. the Lagrangian for spinon is 

L{A, s)~j - A^)s a \ 2 + m 2 s s* a s a j (4.7) 

with /i — 0, 1, 2. It turns out that A^ ~ s* <9 M s a , therefore it is staggered. 

The same gauge field appears in the covariant derivative, obtained via Peierls substi- 
tution, acting on holons. Integrating out spinons and holons one obtains the low-energy 
effective action for A^. Since s a is AF and gapful, the spinon contribution is Maxwellian 
with a thermal mass m for A . On the contrary, holons have a finite FS (in both PG 
and SM) and their contribution to the transverse component exhibits the so-called Reizer 
singularity (Reizer 1989a,b). 

As a result the leading behavior of the A propagator is: with i,j = 1,2 

< AjAj > (u,,q) ~ (S tj - ^(^-Xfy 1 , (4.8) 

where k ~ k F is the Landau damping and x = Xs + Xh with Xs ~ Xh = is the 
total diamagnetic susceptibility, dominated, at least for small S, by the holon contribution. 

There is a characteristic scale of gauge fluctuations emerging from the Reizer singu- 
larity at finite T. It is a kind of anomalous skin depth, derived by assuming as typical 
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energy u ~ T, with the consequence that the transverse gauge interaction is peaked at 



The Reizer singularity in this approach plays a crucial role in the interpretation of 
transport properties, and we remark that it is due to the simultaneous appearance of a 
finite FS and a gap for the bosonic excitations s a ; if these bosons condense it disappears. 
This is what happens at low temperatures in the slave boson approach, which thus has 
the difficulty that the Ioffe-Larkin rule cannot be extended below the holon condensation 
temperature. 

The scalar component Aq has a low energy propagator given by 



In view of the constant term in (14.91) the interaction mediated by Aq is short ranged, 
hence subleading at large distance w.r.t. the interaction mediated by A T , triggered by 
the Reizer (1989a,b) singularity. 



5 The effect of gauge fluctuations on correlation func- 



A difficulty in considering the gauge fluctuations is the lack of a 'small parameter' for 
expansion. Also, a perturbative treatment would be insufficient to get a bound state or a 
resonance with the electron quantum numbers out of a spinon and a holon, whereas some 
kind of binding is expected on the basis of the considerations leading to the Ioffe-Larkin 
rule. A possible way out is to implement the idea of binding using an eikonal approach 
in studying the hydrogen atom as a bound state of a proton and electron. To analyze 
the behavior of correlation functions of physical, hence gauge invariant, fields we apply a 
kind of eikonal resummation of (transverse) gauge fluctuations (Marchetti et al 2004a). 
This resummation is obtained by treating first A^ as an external field, expanding the 
correlation function in terms of first-quantization Feynman paths, then integrating out A^ 
to obtain an interaction between paths which is then treated in the eikonal approximation. 
Finally a Fourier transform is performed to get the retarded correlation function. Further 
approximations are needed however, especially in the treatment of short-scales, to get the 





(4.9) 



tions 
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final result; we refer the reader to Marchetti et al (2004a, 2005) for details, but we briefly 
comment on some interesting features encountered in the calculation of, for example, the 
magnon correlation function. 

5.1 Magnon correlation function 

It turns out numerically that for 5, T small, (identified with PG) the spatial Fourier trans- 
form in the eikonal approximation is dominated by a nontrivial complex \x\- saddle point, 
\x\s.p. ~ QyV 71 "/ 4 due to the effect of gauge fluctuations. The self-consistency require- 
ment for this eikonal approximation yields a region of validity given approximately by 
m sQT ^ T/x < Tn% (inserting numbers as order of magnitude from a few tens to a few 
hundred K). The upper bound temperature, T*, roughly coincides both as order of mag- 
nitude and as 5 — T dependence (for low doping) with the pseudogap temperature, slowly 
decreasing with 5 due to a delicate cancellation, where the relevance of the logarithmic 
correction in (14.61) shows up: T* ~ X m2 s ~ g^|51n<5| ~ ^|ln<5|. The lower bound tem- 
perature T sg , with a maximum at 5 ~ 0.02 and then slowly decreasing, is reminiscent of 
the spin-glass crossover in the heavily underdoped region. 

The main effect of the complex saddle point within the above range is to induce a 
shift in the mass of spinons: m s — ► M = [m 2 s — icT/x) 1 ^ 2 , where c ~ 3 is a constant, 
thus introducing a dissipation proportional to T. Physically, due to the Ioffe-Larkin rule 
electron conductivity is dominated by spinons, and the competition between the mass gap 
and the dissipation appearing in M is responsible for the MIC of in-plane resistivity and 
of many other crossovers in PG, as discussed already in Sec. 2. 

For sufficiently high 5 or T (region identified with SM) one can verify numerically that 
the saddle point contribution is negligible w.r.t. the contribution of fluctuations around 
in the range \x\ < Q^ 1 . The result of gauge fluctuations can also be summarized in 
terms of the appearance of a dissipative term: m s — > m s — idTQT/x m2 j with d ~ 0.1. 
In SM, however, the change of kp from 5 to 1 — 5 yields a decrease of diamagnetic 
susceptibility Xi implying that the thermal de Broglie wave length is shorter than its PG 
counterpart. Therefore the spin-gap effects (£af < ^t) are less effective, being confined 
to very low temperatures. At a pictorial level the existence of a characteristic scale 
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X 



Qt 1 in PG seems to indicate that spinon and antispinon in the magnon have this 
typical size, while in SM QTp 1 is just the typical range of their spatial fluctuations. Our 
approximate time-Fourier transform involves an UV cutoff XQ^ 1 in terms of a parameter 
A fixed phenomenologically by comparison with experimental data. It turns out that 
A << 1 in PG and in SM the result is self-consistent only if rrisQ^ 1 , < 1 with A ~ 1 
(numerically yielding a range from a few tens to a few hundred K ). The retarded magnon 
propagator derived in this manner is given by: 

. . f ^^Jo(|g1(2QT)-V^) PG 

(Q • n>( Wl g) ~ zSM M!i (5.1) 

i2 u>-2m a +i— S-c'Q T A a 

where Jn is the Bessel function, a ~ -^—Qr and <f is measured from the AF wave vector. 

For the wave function renormalization constant we have Z^ G ~ QT^^kpM) 1 / 2 with 
kp ~ <5 and Zq M ~ Q^ 2 kpm s with fcp ~ 1 — 5. Equation (15.11) implies that, within 
the approximation scheme adopted, the transverse gauge fluctuations couple the spinon- 
antispinon pair in Q into a magnon resonance with mass gap 



m n ~ 2m s (5.2) 



and inverse life time 



J SM ~ T if m 2 » ^, ~ T 1 / 2 if m 2 ~ ^ PG 

As briefly discussed below, T^ 1 is proportional to to the electron life-time, which thus 
always increases as T decreases. How then this is compatible with MIC? The point is 
that above MIC the increase of the electron life-time yields an increase of conductivity, but 
below MIC the mobility decreases as T lowers because the spinon can move only through 
thermal duffusion, induced by gauge fluctuations due to gapless holons, and, because of 
binding, the electron can move only if the spinon does. Finally we expect that the effect 
of the neglected residual short-range attraction mediated by Aq gives a renormalization 
of the mass gap, without introducing significative qualitative changes. 
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5.2 Electron Propagator 

Let us now briefly comment on the electron correlation function. Previously we discussed 
the effect of gauge fluctuations on the Q correlation function at large scales, using the first 
quantization path-integral representation for the massive spinon propagator. An analo- 
gous representation is hard to use for the holon correlation function because of the finite 
density of holons. This representation would in fact contain a series of alternating sign 
contributions, corresponding to an arbitrary number of closed fermion wordlines, describ- 
ing the contributions of /i-particles in the finite-density ground state (see e.g. Frohlich and 
Marchetti 1992). To overcome this difficulty, we apply a dimensional reduction by means 
of the tomographic decomposition introduced by Luther (1978) and Haldane (1994). To 
treat the low-energy degrees of freedom we choose a slice of thickness hp /A, with A 3> 1, 
in momentum space around the FS of the holon. We decompose the slice in approximately 
square sectors; each sector corresponds to a quasi-particle field in the sense of Gallavotti- 
Shankar renormalization (Benfatto and Gallvotti 1990, Shankar 1994) (see also Frohlich 
1994). Each sector is characterized by a unit vector n(8), pointing from the center of 
the FS to the center of the box, labelled by the angle 9 between this direction and the 
k x axis. The contribution of each sector can be viewed approximately as arising from a 
quasi ID chiral fermion; this avoids the finite FS problem for the path-integral represen- 
tation discussed above. The paths appearing for a sector are straight lines directed along 
the Fermi momenta of the sector, with small gaussian transverse fluctuations. We now 
multiply this path representation for the holon by that for the spinon, following the de- 
composition formula (13.81) for the electron treated in MF. Finally we integrate out in the 
eikonal approximation the transverse gauge fluctuations, proceeding as for the magnon 
correlation function. The output for the retarded correlation function becomes simple 
only for momenta on the FS, however the main features presented here hold in a slab of 
scale Qt around the FS. On the FS the result can be summarized as follows: 



In this approximation the scattering rate T comes entirely from the spinon, hence it 
is half of Tq. The wave-function renormalization constant comes as a product of a spinon 




(5.4) 



19 



and a holon contribution plus the effect of the quasi ID structure of holons yielding for 
(the isotropic component of) Z: 

Z ~ (Qrm,) 1 / 2 . (5.5) 

By chiral invariance, physical absence of superconducting gap, the mass term induced by 
the spinon should be exactly cancelled by the holon-spinon short range attraction. Notice 
that the same assumption underlines the recovery of the FS in the SU(2) slave-boson 
approach (Lee et al 2006). We should emphasize that the electron FS , coinciding in our 
treatment with the holon FS, in PG is the 'small FS' centered at (±f , ±f ) with fcp ~ 5 
whereas in SM it is the 'large FS' with kp ~ 1 — 5. 

In the PG phase besides (15. 5p we have an additional angular dependent factor around 
e.g. (~, ~) given by — ^[cos(arg /cp)+sin(arg k F )}} of the wave function renormalization 
yielding a reduction of the spectral weight outside the magnetic Brillouin zone, reminiscent 
of the 'Fermi arcs' of ARPES in underdoped cuprates ( see Wen and Lee (1996), Lee 
et al (1998) for a similar situation in the SU(2) x U(l) slave boson approach). This 
anisotropy is due to the matrix structure of the 'Dirac' hamiltonian, describing physically 
the 'zitterbewegung', and the relation between the original spinless holon and the two 
species of 'Dirac holons' introduced by Bmf-, which changes in the four quarters of the 
Brillouin zone centered at (±f, ±f ), thus recovering the standard translational invariance 
for the electron. We conjecture (Marchetti et al 2007) that an additional reduction of 
the spectral weight appears as a consequence of the attraction between spin vortices of 
opposite vorticity, neglected in the present MF, improving the agreement with ARPES 
data. 

The structure exhibited by (15. 4p shows that the gauge fluctuations are able to bind 
together spinon and holon into a resonance for low energies at the Fermi momenta, but 
with a strongly temperature dependent wave function renormalization constant. In par- 
ticular Z ~ T 1//6 , so Z vanishes if formally extrapolated to T = 0. This implies a peculiar 
non- Fermi liquid character for this system of 'electron resonances'. In fact, by gauge 
invariance the sum of spinon and holon currents is zero, but at T = the spinons can- 
not move as they are gapped and without any possibility of thermal diffusion, therefore 
the electron system behaves as an insulator. A similar phenomenon is discussed within 
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Weng's approach in Kou and Weng (2005). Existence of a FS is then compatible only 
because Z = 0. However, a real extrapolation to T = cannot be done because of the 
range of validity of the approximation, as discussed above. The system therefore appears 
to fit naturally within the scheme of unstable fixed points (UFP) outlined by Anderson 
(2002). 

6 Transport properties 

In this section we sketch the results of calculation for physical transport quantities through 
Kubo formulas using the above approach and compare the theoretical results with the 
experimental data. 

6.1 In-plane resistivity 

The in-plane resistivity p is calculated via the Ioffe and Larkin (1989) addition rule: 
p = p s + ph, where p s is the resistivity of the spinon-gauge subsystem and pt of the 
holon-gauge subsystem, which is subdominant and of standard form for a gauge-fermion 
system (Nagaosa and Lee 1990, Lee and Nagaosa 1992). Via the Kubo formula: 

(p,)- 1 =2f Q oo dx°x (ff)(x ,q = 0), (6.1) 

where j s ~ dQ ~ Qt& is the spinon current and, using (15.11) one obtains in PG: 

From (16. 2p one recovers the MIC and, due to the square root in M , an inflection point 
T* ~ x m2 s a t higher temperature, found also experimentally (Takagi et al 1992). Hence, in 
our approach the MIC is due to correlation effects, not to a disorder-induced localization, 
as discussed in Sec. 2. In fact, within the range of validity for deriving ( 16.21) (above the 
Neel temperature and the 'spin-glass' temperature) the insulating behavior is power-law 
like, not exponential. We believe that outside this range the qualitative behavior should 
be the same, but the formula itself cannot be applied directly. 

Moreover, the normalized resistivity p n = (p — p(Tmic))/{p(T*) — p(Tmic)) is a func- 
tion only of the ratio y = cT/T* : p n (y) — (1 + y 2 ) 1 ^ 8 / sin(| arctany). A universal 
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behavior of this quantity in fact can be verified in experimental data and in general terms 
has already been empirically noted in publications (Wuyts et al 1996, Konstantinovic et al 
2000). Furthermore, the dependence of \M\ on T 2 yields an approximate linear behavior 
of in-plane resistivity in T 2 , for temperatures just above the minimum in PG (Marchetti 
et al 2006), observed experimentally (Ando et al 2004), without involving 'Landau' quasi- 
particles advocated in (Ando et al 2004). 

The introduction of a magnetic field perpendicular to the plane allows one to find MIC 
also for higher dopings suppressing superconductivity, and it yields a shift of the MIC to 
higher temperature which in turn produces a big positive transverse magnetoresistance 
(Marchetti et al 2001), as found in experiments in LSCO (Lacerda et al 1994, Kimura et 
al 1996, Abe et al 1999). 

Non-magnetic impurities simply affect the resistivity by an upward shift as they ap- 
pear via the Matthiesen rule only in ph, consistent with irradiation experiments (Rullier- 
Albenque et al 2003). In contrast, Zn doping provides strong potential scattering changing 
the holon resistivity, and at the same time disturbs the AF background, making the AF 
correlation shorter, therefore qualitatively we expect a shift of MIC in spinon resistivity 
up in T, as observed in Zn-doped YBCO (Segawa and Ando 1999). This strong contrast 
of behavior for non-magnetic (affecting the holon contribution to resistivity) and magnetic 
(affecting the spinon contribution) scattering, is indirect proof for the Ioffe-Larkin rule: 
It is the resistivity, not the conductivity of the constituent slave particles which is added 
to form the resistivity of the physical system. 

On the other hand, in SM one obtains : 

p ,,2AQ f '(r + ^/r) = ^ + ig, (6.3) 

where A ~ 0.7. In the high temperature limit Qt 3> m s , the damping rate in ( 16.31) 
dominates over the spin gap 2m s and the spinon contribution to resistivity is linear in 
T, with a slope a ~ (1 — 5)/(<5| ln<5|). Lowering the temperature, the second term in 
(16.31) gives rise first to a superlinear behavior and then, at the margin of validity of our 
approach, an unphysical upturn. The deviation from linearity is due to the spin gap 
effects and is cut off in the underdoped samples by the crossover to the PG phase. We 
expect that physically in the over doped samples it is cut off by a crossover to a FL' phase'. 
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The celebrated (approximate) 'T-linearity' can also be understood qualitatively as 
a consequence of T and the 'effectiveness' of the gauge fluctuations to form resonance, 
predominantly contributing to conductivity (see previous Section), in a slab of momenta 
Qt ~ T 1 / 3 around the Fermi surface. In fact the conductivity derived from the Boltzmann 
transport theory would be o~o ~ T -1 ) but due to 'effectiveness' the physical conductivity 
is a ~ cfqQt ~ T _4 / 3 T 1 / 3 ~ T _1 . Similar considerations apply to the linearity found in 
SM for out-of-plane resistivity, spin relaxation time 63 (TTi) and, replacing T by uo of the 
AC conductivity (Marchetti et al 2005), as discussed below. 

6.2 In-plane IR electronic AC conductivity 

We compute the electronic AC conductivity in the two limits u << T, T << uj, where 
uo is the external frequency. It turns out that up to the logarithmic accuracy one can 
pass from the first to the second limit by replacing T with uj in Qt and T ( we denote 
the obtained quantities by Q w , TJ) and rescaling T by a positive multiplicative factor 
A < 1/2. In PG, in the range T < uj < Jm s the AC conductivity can be approximately 
obtained from p _1 derived from eq. fl 6 . 2 H substituting T with Xu. Hence it exhibits a 
broad maximum at a frequency ujmic (self-consistently << Jm s ), corresponding to a 
temperature slightly higher than Tmic, as in fact experimentally seen in Dumm et al 
(2003). An a — b anisotropy of T MIC and ujmic found in underdoped LSCO (Dumm et 
al 2003) would be generated naturally in the above scheme by an a — b anisotropy of the 
magnetic correlation length (Marchetti et al 2004b). The anisotropy of Tmic strengthens 
our interpretation of correlation- induced MIC, since a disorder-induced localization is 
expected to have a unique MIC temperature. As T becomes greater than u the AC 
conductivity becomes approximately uj independent. 
In SM, for uo « T we have 



a(u,T) 



Qt 1 



(6.4) 



i(uj - 2m s ) + T i(u- 2m s )T- 1 / 3 + T ' 



while for T << uj 



a(u,T) 



CI 



1 



(6.5) 



i(uo — 2m s ) + Ar, 



i(uj — 2m s )w -1 / 3 + uj 
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From the above formulae the following features found experimentally in overdoped LSCO 
(Startseva et al 1999) and BSCO (Lupi et al 2000) are easily derived: besides the standard 
tail ~ u' 1 , the effect of replacing Qt by Q u is to asymmetrize the peak at 2m s , appearing 
in Kea(u,T) for uj « T and to shift it towards lower frequency. 

6.3 Spin-lattice relaxation rate 63 (l/XiT) 

Assuming Qt as the cutoff for the | q\ integration and using the smoothness of the hyperfme 
field at this scale we derive in PG: 

{a + bT m 2 » — 
s X 

with a/6 ~ 0.1, and the two terms are due to Re and Im of J in (15.11) . 

Therefore we obtain a broad peak as observed in some cuprates (Berthier et al 1994) 
and up to a multiplicative constant a universality curve as a funcion of y = cT/T* can be 
derived. 

In SM one finds 

8S (j7y)~(l-*) 2 P7 1 - (6-6) 
Therefore, in the high temperature limit we recover the linear in T behavior for 63 (T\T) 
and at higher dopings /lower temperatures the superlinear deviation, also found experi- 
mentally in overdoped samples of LSCO (Berthier et al 1997, Fujiyama et al 1997). 

6.4 Out-of-plane resistivity 

We assume an incoherent transport along c-direction. p c thus is dominated by virtual 
hopping between layers and it is calculated via the Kumar- Jayannavar formula (Kumar 
and Jayannavar 1992, Kumar et al 1997, Kumar et al 1998): 

where t c in our approximation is an effective average hopping in the c-direction and v the 
density of states at the Fermi energy. In PG the first term in ( 16. 71) dominates, yielding 
an insulating behavior: 



p c ~ (|M| sin(argM))" 



T 2 ml 



,2 cT 
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This crossover reproduces a knee experimentally found in Ito et al (1991) and Yan et al 
(1995). A very strong decrease with T of the anisotropy ratio p c {T) / p ah {T) in PG is also 
well reproduced by the above formulae. 

In the SM phase the second, metallic, term in (16.71) dominates and substituting V and 
Z one recovers the T— linearity in the 'incoherent regime' V ^> t c Z\ 

j2 T ~ T 
Pc ~ t 2 c m s u X m 2 s ~ (S\lnS\y/ 2 ' (6 ' 8) 

A universal behavior is also recovered for underdoped samples in the form p c = a c (5)p c (y) 
with p c (y) ~ (1 + y)~ 1 ^ j sin(| arctan?/) in PG, ~ y in SM and the S dependence of a c is 
in agreement with the phenomenological proposal of universality (Su et al 2006). 



7 Concluding Remarks 

As mentioned in the Introduction, there is no consensus theory of HTS cuprates yet. What 
has been briefly summarized in this paper is a theoretical attempt to describe the metal- 
insulator crossover as the most spectacular phenomenon of cuprates in the pseudogap 
phase, as well as other transport properties of cuprates in the pseudogap and 'strange 
metal' phases in a self-consistent way, using the spin-charge gauge approach. To the best 
of our knowledge, this is the only available theoretical treatment of this crossover as an 
intrinsic property of the doped Mott insulators. We start from a two-dimensional t — J 
model in the limit t > J to describe the strong correlation effects and use the gauge 
field approach to implement the single-occupancy constraints. In the mean field theory 
we have chosen, the spin degrees of freedom are represented by hard core bosons, while 
the charge degrees of freedom are represented by spinless fermions, with a slave particle 
gauge field binding them 'loosely'. 

We have compared the calculated results with experiments, and in comparing the 
temperature-doping dependence of the in-plane resistivity, including the metal-insulator 
crossover (Marchetti et al, 2000, 2001, 2004a) there are no adjustable parameters, ex- 
cept for the absolute scale of resistance. Moreover, the nontrivial doping/magnetic field 
dependence of the metal-insulator crossover temperature is a direct consequence of the 
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theoretically predicted spinon mass gap m 2 s ~ <5j log <Jj , where 5 is the doping concen- 
tration, and the competition of the short-range antiferromagnetic order with dissipative 
motion of the charge carriers. Using the metal-insulator crossover as a clue to understand- 
ing the nature of the pseudogap phase, we could explain a number of peculiar transport 
properties in the pseudogap phase including the 'universality' and anisotropy of in-plane 
conductivity, the huge positive magnetoresistance in underdoped samples, the 'knee' in 
c-axis resistivity, the broad peak in the nuclear relaxation rate and low-frequency spin 
fluctuation spectrum, within the same framework, using the same set of parameters. 

Using the 'melting' of the n flux lattice as the main signature of crossover from the 
pseudogap to the strange metal phase, we could easily generalize our approach to higher 
doping/temperature and calculate various transport properties to compare with exper- 
imental data in this phase. Apart from recovering the well-known linear temperature 
dependence of in-plane and out-of-plane resistivity, NMR relaxation time, u' 1 asymp- 
totic behavior of the AC conductivity, etc, we could also account for some more subtle 
effects like the deviation from linearity and the asymmetry of the AC conductivity peak 
(Marchetti et al 2005). 

Within the gauge field approach it is natural to consider the low-energy excitations 
as 'composite particles' consisting of slave particles loosely bound by the gauge field. A 
question arises: Is this gauge field a pure mathematical construction or does it reflect to 
some extent the physical reality? In our opinion, the second is true. In strongly correlated 
systems the spin and charge degrees of freedom are neither 'completely separated' nor 
'bound in a point-like particle', so the gauge field is an 'indispensable companion' of 
slave particles describing the charge and spin degrees of freedom. It took some time for 
condensed matter physicists to accept the quasi-particle as a real physical entity. The 
quasi-particle is a constituent of collective motion, and it ceases to exist as soon as it is 
taken away from the condensed matter. In a similar way, the gauge field disappears as 
soon as the strong correlation is switched off. In this sense, it is reasonable to consider 
low-energy excitations in strongly correlated systems as composite particles bound by the 
gauge field. 

We would like to sincerely thank J.H. Dai, L. De Leo, G. Orso and F. Ye for their 
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collaboration at different stages of this project. 
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